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Bond-dilution effects on the ground state of the square-lattice antiferromagnetic Heisenberg 
model, consisting of coupled bond-alternating chains, are investigated by means of the quantum 
Monte Carlo simulation. It is found that, when the ground state of the non-diluted system is a 
non-magnetic state with a finite spin gap, a sufficiently weak bond dilution induces a disordered 
state with a mid gap in the original spin gap, and under a further stronger bond dilution an 
antiferromagnetic long-range order emerges. While the site-dilution-induced long-range order is 
induced by an infinitesimal concentration of dilution, there exists a finite critical concentration 
in the case of bond dilution. We argue that this essential difference is due to the occurrence 
of two types of effective interactions between induced magnetic moments in the case of bond 
dilution, and that the antiferromagnetic long-range-ordered phase does not appear until the 
magnitudes of the two interactions become comparable. 

KEYWORDS: dilution-induced antiferromagnetic long-range order, quantum Monte Carlo simulation, 
quantum phase transition, non-magnetic impurity, spin-Peierls compound 



1. Introduction 

Randomness effects on quantum antiferromagnetic 
(AF) Heisenberg models with a spin-gapped ground 
state have attracted much interest in relation to the 
impurity-induced AF long-range order (LRO) observed 
experimentally, e.g., in the first inorganic spin-Peierls 
compound CuGeOs, 1 ' 2 the two-leg ladder compound 
SrCu2C>3, 3 and the Haldane compound PbN^A^Os- 4 In 
the spin-Peierls compound CuGeC>3, the lattice dimcr- 
ized state associated with a formation of spin gap is 
realized at low temperatures. 1 When a small amount 
of non-magnetic impurities Zn or Mg are substituted 
for Cu, the spin-gapped ground state of the pure sys- 
tem changes to an AF LRO state, 5-7 thereby the lat- 
tice dimerization is preserved. Such an impurity-induced 
AF LRO has been observed also in the bond-disorder 
system, CuGei-^Si^Oa. 2,8 ' 9 In early theoretical and nu- 
merical studies on quasi-one-dimensional diluted Heisen- 
berg antiferromagnets, the effect of inter-chain interac- 
tion is often taken into account by making use of mean- 
field approximations, where the site and bond dilutions 
are unable to be distinguished with each other. 10-12 By 
the recent numerical analyses, however, the characteris- 
tic feature of the site-dilution-induced AF LRO has been 
understood more clearly by treating the inter- and intra- 
chain interactions on an equal footing. 13-15 More recently 
we have carried out extensive numerical simulations also 
on the bond-diluted system, and have found that the 
mechanism for the AF LRO to appear in this system is 
essentially different from that in the site-diluted system. 
The purpose of the present paper is to discuss our results 
in more detail, which have been reported in refs. 16 and 
17 briefly. 

In the present work we concentrate on the two- 
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Fig. 1. Ground state phase diagram of the non-diluted 2D S = 
1/2 AF Heisenberg model. The circles denote the phase boundary 
between the spin-gapped and AF LRO phases. The data points 
except for a = 0.5 were estimated by the previous quantum 
Monte Carlo simulation. 20 The error bar of each point is much 
smaller than the symbol size. The dotted and dashed lines denote 
a = 0.5 and J' = a, respectively. 



dimensional (2D) S =1/2 AF Heisenberg model consist- 
ing of coupled bond-alternating chains. The magnitude 
of the stronger (weaker) intra-chain interaction is put 
unity (a) and that of the inter-chain interaction J' (see 
also Fig. 2 below). The ground state of decoupled chains, 
i.e., J' = 0, is the dimer state with a finite spin gap, A p , 
except at the uniform point (a = l), 18 which has a criti- 
cal ground state with the Gaussian universality. 19 As J' 
increases, the quantum fluctuations, which invoked the 
dimer state, are gradually suppressed and the spin gap 
decreases. When J' exceeds a certain critical value, J' cv , 
which is a function of a, the spin gap vanishes and the 
AF LRO emerges. A global ground state phase diagram, 
obtained by the previous quantum Monte Carlo Simula- 
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Fig. 2. Effective spins (arrows) and an effective AF interaction in- 
bctwccn (Jmn) in a site-diluted system. The thick and thin lines 
denote the strong and weak interactions, respectively, and the 
small circles represent the removed sites. The circles surrounding 
the arrows describe the extent of the effective spin schematically, 
which radius corresponds to the correlation length of the pure 
(i.e., non-diluted) system. 



tion, 20 is presented in Fig. 1. 

The ground state in the spin-gapped phase is described 
qualitatively well by a direct product state of singlet 
dimers sitting on each strong bond, which is called the 
valence-bond solid state. When spins are randomly re- 
moved (site dilution) , a spin which formed a singlet pair 
formerly with a removed spin becomes nearly free. We 
call it an effectively free spin, or simply an effective spin. 
Due to the presence of weak interactions with the sur- 
rounding spins, an effective spin is somewhat blurred] 
its linear extent is proportional to the correlation length 
of the non-diluted system. Between two effective spins 
at sites m and n, there exists an effective interaction, 
J mm mediated by the sea of singlet pairs as illustrated 
in Fig. 2. It is cither ferromagnetic or AF depending if 
m and n are on the same sublattice or on the different 
sublattices. It preserves its staggered nature with regard 
to the original square lattice and thus the dilution does 
not introduce frustration. In this sense we call {J mn } 
the effective AF interactions between effective spins. Its 
magnitude is finite, though exponentially small, even if 
x, the concentration of site dilution, is vanishingly small. 
Consequently, we expect x^., the critical concentration of 
site dilution above which the disorder-induced AF LRO 
appears, is strictly zero. 13-15 

When a strong bond of amplitude unity is removed, 
on the other hand, two spins on the both edges of the 
removed bond become nearly free as shown in Fig. 3. 
We can also notice from the figure that there exists an 
effective AF interaction, J a f, between these two spins, 
which is mediated mainly through the shortest paths of 
the interactions connecting them, i.e., J a f ~ J' 2 . This in- 
teraction alone recombines the two spins to form a singlet 
pair with an excitation gap A, which is proportional to 
J a f. It yields a localized low- lying excited state, so-called 
mid- gap state, as long as A <C A p . In the case of random 
bond dilution, in addition to J a f , there also exist effective 
AF interactions, { J mn \, between effective spins on edges 
of different removed bonds as in the site-diluted system 
(Fig. 3). Naturally, we expect competition between the 



Fig. 3. Two types of effective interactions J m n and J a f in the 
bond-diluted system. The former are mediated by the sea of sin- 
glet pairs as the same as the site-dilution case, whereas the latter 
act through the shortest 2D paths. 



two effective interactions, J a f and { J mn \- 

In the present paper, we discuss how the competition 
between these two types of interactions gives rise to the 
bond-dilution-induced AF LRO associated with a finite 
critical concentration of bond dilution, x c (> 0), in con- 
trast to x s c = for the site dilution. We also investigate 
in some detail the mid-gap state due to J a f in the weak 
dilution limit, x — > 0, and shortly refer to a peculiar 
phenomenon which originates from the randomness in 
the system, namely, a possible existence of the quantum 
Griffiths phase neighboring to the AF LRO state. 

The present paper is organized as follows. In §2, the 
model and the method of our numerical analyses are in- 
troduced. Detailed analyses of the mid-gap state in the 
limit of x — > are presented in §3, and in §4 the ground 
state phase diagram on the x-J 12 plane is determined. 
Section 5 is devoted to summary and discussion. 

2. Model and Numerical Method 

We consider the bond-diluted quantum AF Heisenberg 
model on the square lattice of coupled bond-alternating 
chains. Its Hamiltonian is written as 

W = /, e (2i,j)(2i+l,j)^2i,j ■ Sai+lJ 



a e (2i+l,j)(2i+2,j)^2i+l,j 1 ^2i+2,j 



hJ 



H,j+1 > 



(1) 



where 1 and a (> 0) are the AF intra-chain alternat- 
ing coupling constants, J' (> 0) the AF inter-chain cou- 
pling constant, and Sij is the quantum spin operator 
with magnitude S = 1/2 at site We choose the 

x-axis as the one along chains and the y-axis as in the 
inter-chain direction. Randomly quenched bond occupa- 
tion factors {t(i,j)(k,l)} independently take either 1 or 
with probability 1 — x and x, respectively, where x is the 
concentration of bond dilution. 

The quantum Monte Carlo (QMC) simulations with 
the continuous-imaginary-time loop algorithm 21 ' 22 are 
carried out on L x L square lattices with periodic bound- 
ary conditions. For each sample with a bond-diluted con- 
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Fig. 4. J' dependences of the the correlation lengths and the 
spin gap A p in the non-diluted system with a = 0.5, where 
Jc,p = 0.5577(1) is the critical value of J' determined by the 
same method as in ref. 20. The dotted line denotes J' = 0, the 
system of decoupled chains. The error bar of each point is smaller 
than the symbol size. The solid and dashed lines going through 
the data points are guides to eyes. 
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Fig. 5. Dependences of A' 2 ' and A( 4 ) estimated by the second- 
and fourth-moment estimators on x = 1/(2L 2 ) for a = 0.5 and 
L X L systems with L = 16, 18, 20, 24, 32, 48, 64, and 80. The 
lines are the fitting functions [eqs. (5) and (7)]. 



figuration, 10 3 ~ 10 4 Monte Carlo steps (MCS) are spent 
for measurement after 5 x 10 2 ~ lfj 3 MCS for thermal- 
ization. For each x the random average is taken over 
10 1 ~ 10 3 samples. As the temperature is decreased, each 
measured quantity, e.g., the structure factor, converges 
to a finite value. This reflects a finite excitation gap due 
to either of an intrinsic spin gap or the finiteness of the 
system. The simulations are performed at low enough 
temperatures, typically T ~ 10~ 3 , where all the physical 
quantities of our interest do not show any temperature 
dependence and thus can be identified with those at the 
ground state. 

Before going into discussions on bond-diluted systems, 
let us here mention some ground state properties of the 
non-diluted systems described by Hamiltonian (1) with 
x = 0. Its ground state phase diagram on the a- J' plane 
is shown in Fig. 1. For discussions below we show in Fig. 4 
the J' dependence of the spin gap A p and the inverse ge- 
ometric mean of the correlation lengths £p and ^ along 
a = 0.5, on which we concentrate our analyses of bond- 
dilution effects below. The correlation lengths £p and 
£v are estimated from the dynamical correlation func- 
tions at momenta (w, 0) and (0, w), respectively, by using 
the second-moment method. 23 Similarly, A p is evaluated 
from the correlation function along the imaginary-time 
axis (see also below). As seen in Fig. 4, both A" 1 and 
£p£ p / exhibit power law divergence as J' approaches J' c 
from below. 

3. Mid-Gap State in the Limit of x — * 

When a strong bond is removed from the non-diluted 
system with the non-magnetic ground state, two mag- 
netic moments induced at the both ends of the diluted 
bond are expected to reform a spin singlet through J a f, 
one of the effective AF interactions mentioned in §1. For 
a sufficiently weak bond dilution, x <^ 1, we expect a 
mid-gap state to appear, since the spin gap A induced 
by J a f is smaller than A p , the spin gap in the non-diluted 



system. In this section, based on the QMC data, we dis- 
cuss the J' dependence of the gap attributed to J a f in 
the limit of x — > 0. 

To extract the value of A, we make use of the second- 
moment formula, 22 ' 23 



' C(0) 
<5(2tt//3) 



(2) 



Here (3 is the inverse of temperature and C(u>) — 
Jq drC(r)e 4WT is the Fourier transform of the imaginary- 
time staggered correlation function 

c ( T ) EE Z2E(" 1 ) |r, " rjl ^(°)^( r )) ■ 



(3) 



In the bond-diluted system of present interest, it is ex- 
pected that C(t) is dominated by two terms; the one 
representing the original spin gap in the non-diluted sys- 
tem and the one attributed to the mid gap. That is, the 
imaginary-time staggered correlation function can be ap- 
proximated as 17,24 



C(r) 



cosh 



cosh 



2 
0A 



cosh{(r-/3/2)A p } 
1 cosh{(r-/3/2)A} , (4) 



where a p and Oj are some functions of a, J', and x. The 
coefficient a u which represents the weight of the contri- 
bution from the mid-gap state, is naturally considered 
to be proportional to the concentration of diluted bonds, 
i.e., Oi = ax with a = a(a 7 J'). From eqs. (2), (4), and 
a D + a; = 1, A^ 2 ) is evaluated as 



A< 2 > = AJ(1 



1 — ax A 



)/{l 



1 



II UiiXj L — ip Lilly . — ip 

for <C A <C A p . Note that in the last expression, 
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Fig. 6. J' 2 dependences of the spin gap estimated by the second- 
and fourth-moment estimators. The dashed line, which repre- 
sents A = 1.5J' 2 , is a guide to eyes. 



A^ 2 ) explicitly depends on x as well as on A. We regard 
eq. (5) as a fitting function to extract A in the weak 
dilution limit, where A is assumed to be independent of 
x. 

The QMC simulation to evaluate A^ 2 ) is performed 
along a = 0.5 for the system with only one diluted bond, 
which corresponds to the x = 1/2L 2 system. By taking 
the limit of L — > oo, we can investigate the properties in 
the limit of x — > 0. The obtained set of A( 2 )(a;) for each 
J' is fitted by eq. (5), thereby the value of A p is fixed to 
the one simulated in the non-diluted system. In the upper 
part of Fig. 5 the x dependence of A^ 2 ) is presented. The 
resultant spin gap A is shown in Fig. 6. As a function of 
J', it increases proportionally to J' 2 , while A p decreases 
(see Fig. 4). 

We have further confirmed the values of A by using 
the fourth-moment estimator, 22 defined by 



(A«)- 




C(0) - C(2w/0) 



C(2tt//3) - C(4tt//J) 
In this case the fitting function becomes 



1 



(6) 



aw =aJ{i 



1 



^(^) 3 }/{l 

ax A n 



1 



ax A n 



(7) 

The lower part of Fig. 5 is the x dependence of A^ 4 ) 
obtained by eq. (6). It is observed that as x increases 
A^ 4 ) converges more rapidly to a constant values, which 
corresponds to A. This is due to the smaller corrections, 
[(1 - ax)/ax]{A/A p ) 3 1 in eq. (7) than those for A< 2 ). 

Figure 6 shows a nice coincidence of the values of A, es- 
timated by the second- and fourth-moment methods. We 
therefore conclude the existence of a mid-gap state whose 
gap A is proportional to J' 2 and distinctly smaller than 
A p , and whose contribution to the imaginary-time cor- 
relation function is proportional to x. This is just what 
we have expected, i.e., a reformed singlet pair due to the 
effective AF interaction J a f acting on two effective spins 
at both ends of a diluted strong bond. 



Fig. 7. System size dependences of S S (L, 0, x)/L 2 for x = 0.02 
(crosses), 0.04 (x-marks), 0.06 (squares), 0.08 (circles), 0.1 (up- 
ward triangles), 0.2 (downward triangles), 0.3 (diamonds), and 
0.35 (pentagons) at (a, J') = (0.5,0.5). Dashed lines are ob- 
tained by the least-squares fitting. The extrapolated values are 
denoted by the filled symbols which correspond to A£? /3 in the 
thermodynamic limit. 
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Fig. 8. Dependences of the staggered magnetization on the con- 
centration of bond dilution for (a, J') = (0.5,0.5) and (0.5,0.3). 
The solid circle and triangle denote the quantum critical points, 
x c ss 0.0502 and 0.1101, respectively, estimated by the finite- 
size scaling analysis described in the text. For comparison, we 
also show the result of the site dilution for (a, J') = (0.5,0.5), 
which is taken from ref. 15. The percolation thresholds on the 
site and bond processes are denoted by x^ ss 0.41 and ijj = 0.5, 
respectively. All the lines are guides to eyes. 



4. Bond-Dilution-Induced Antiferromagnetic 
Long-Range-Ordered Phase 

As x is increased from zero, fitting of C(t) by us- 
ing eq. (4) becomes unstable, implying that a number 
of mid-gap states, other than A discussed in the previ- 
ous section, appear. We have clearly observed, however, 
that A^ 2 ) of eq. (2), which we now regard as an upper 
bound of the spin gap, as well as the inverse of the spatial 
correlation lengths, £ x and £ v , decrease as x increases. 17 
Finally, when x exceeds a certain critical value, x c , the 
AF LRO appears. During the present section, we call 
the state between a; = and x c simply the disordered 
state, and concentrate on the quantum phase transition 
at x = x c to the AF LRO state, postponing the discus- 
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Fig. 9. Finite-size scaling plots of the correlation length for (a) (a, J') = (0.5,0.5), (b) (0.5,0.4), (c) (0.5,0.3), and (d) (0.5,0.2). 



sion about a possible existence of the quantum Griffiths 
state at x < x c until the last section. 

To discuss the quantum phase transition between the 
disordered and the AF LRO phases, first we investigate 
the zero-temperature staggered magnetization, M s , as a 
function of x, which is evaluated as 

M*(x) = lim lim 3S °( L ' :7 » . ( 8 ) 

Here S s (L,T,x), the static structure factor at momen- 
tum (jt, 7r), is defined by 

S s (L,T,x) = ±.J2(-i)\r*-r>\(SfS3) . (9) 

In Fig. 7 we show the system size dependence of 
S a (L,0,x)/L 2 for various x's. For x > 0.06, the QMC 
data for each x are fitted fairly well by a linear expres- 
sion, a + bL~ x , where the value of a, represented by the 
solid symbols in the figure, gives an estimate for M s 2 /3 in 
the thermodynamic limit. For x = 0.02 and 0.04, on the 
other hand, we find a linear fit extrapolates to a negative 
il/ s 2 /3 for L — > oo, which indicates a vanishing staggered 
magnetization. Indeed, the data are fitted much better by 
S S (L, 0, x) = a(l — cxp(— bL)), which is the scaling form 
derived for a non-magnetic ground state in the modified 
spin- wave theory. 25 

In Fig. 8 the x dependences of M s arc shown for 
J' = 0.3 and 0.5, together with that for the site-diluted 
system 15 for comparison. The value of a is fixed to be 



0.5 in all the cases. It is clearly demonstrated that while 
the AF LRO is induced at an infinitesimal concentration 
of site dilution, in the bond-diluted system there exists 
a critical concentration, e.g., x c w 0.11 for J' = 0.3. It is 
also definitely finite even for J' = 0.5, whose non-diluted 
system is close to the critical point J' c m 0.5577. 

In passing we note that M s (x) shown in Fig. 8 has a 
peak at x = x opt . In the strong dilution regime, x op t ^ 
x < Xp' b , the further dilution tends to destroy the AF 
LRO and the latter would vanish just at the percolation 
threshold. 26 ' 27 

Next let us determine the phase diagram on the x-J' 2 
plane. The critical concentration of dilution, x c , which is 
roughly evaluated above, is estimated more precisely by 
the finite-size scaling analysis of the correlation length 
£ x as shown in Fig. 9. The critical concentrations are 
estimated as x c = 0.0502(2), 0.0926(2), 0.1101(2), and 
0.1161(6) for J' = 0.5, 0.4, 0.3, and 0.2, respectively. 
[Almost the equal values of x c are obtained by the same 
analysis on £ v .) The values of the critical exponent of the 
correlation length are v = 0.69(5), 0.74(2), 0.84(3), and 
0.95(4), respectively, which monotonically increases from 
v = 0.71, the critical exponent of the three-dimensional 
classical Heisenberg universality class. 28 We note, how- 
ever, that the finite-size scaling analysis with the fixed 
critical exponent v = 0.71 gives the same value of x c 
within the numerical accuracy. Certainly, more precise 
numerical data are needed to settle this point. 

The resultant phase diagram is shown in Fig. 10. In 
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x 

Fig. 10. Phase diagram of the ground state for a = 0.5. The cir- 
cles are evaluated by the finite-size scaling analysis. The dashed 
line is a guide to eyes. In the inset we show the relation be- 
tween l/y^CpSp^'c and In J' 2 , where we used the QMC data 

(£p.£p) = (11-998(9), 9.312(10)), (5.320(1), 3.460(2)), (3.452(1), 
1.8202(3)), and (2.5492(4), 1.03244(5)) for J' = 0.5, 0.4, 0.3, and 
0.2, respectively. The straight line in the inset represents a lin- 
ear fitting of the data by y = a + bx with a = —0.633(6) and 
b = -0.762(3). 

the inset of the figure we show an important result, i.e., 
the relation of x c to J' 2 , or the inverse of the function 
J' c (x), a critical value of J' for a given x (and a given a). 
Note that J' cp = J' c {x = 0). The effective AF interaction, 
J mn , introduced in §1 is approximately described by 

Jmn oc (-l)l r ">- r " +1 lexp( jL=) , (10) 

where I = \r m — r n \ is the distance between the effec- 
tive spins 29,30 at the ends of different removed bonds. 
On the other hand, the effective AF interaction, J a f , be- 
tween effective spins on the same removed bond is simply 
described by 

Jaf J' 2 ■ (11) 

Thus the relation indicated in the inset of Fig. 10 is read 
as 

J{X C ) Jaf , (12) 

where J(x) is an average of {J mn } acting on effective 
spins on neighboring diluted bonds, and is approximately 
represented as 

,7(3?) ~ Jo exp( t==) > (13) 

in which the exponential part is obtained from J mn of 
eq. (10) with 1/ ^Jx, i.e., the mean distance of the nearest 
removed bonds, substituted for I. The result of eq. (12) 
implies that the phase transition between the disordered 
and AF LRO phases occurs when the magnitudes of Jaf 
and J{x) become comparable. We interpret this result 
as described below. 

Let us replace a bond-diluted system of present inter- 
est by a model system consisting only of effective spins 
with interactions {J m „} and J a f, neglecting all spins in 



the sea of singlet pairs. Since {J mn } do not introduce 
frustration at all, combined with the fact that { J mn } are 
short-ranged except for the system with J' = J^ p , the 
actual geometry of { J mn } is not relevant to the present 
problem, and so we may further simplify the model sys- 
tem to a regular square array of the effective spins; along 
the chain direction the interactions are alternating be- 
tween J a f and J7"(x), while in the inter-chain direction 
the interaction is uniform, i.e., also J{x) given above. 
In this way, the original system in the disordered state 
can be approximately mapped to a gapped state of the 
original non-diluted system of Fig. 1 but with both a 
and J' replaced by J{x)/J a f. As x (and/or J 1 ) of the 
bond-diluted system increases, 3{x)j J & i docs so, which 
corresponds to the upward movement of the mapped reg- 
ular system along the dashed line (J' — a) in the figure. 
Finally it hits the critical line at x = x c (J') (J' = J' c (x)), 
where eq. (12) holds, and the quantum phase transition 
to the AF LRO phase occurs. 

Our above scenario can explain the behavior of the 
critical line in the x — * limit, i.e., J' c {x) at least first 
decreases from J' when the dilution is introduced (see 
Fig. 10). Naively one might expect the dilution sup- 
presses the bulk correlation represented by £p' y , and in 
order to recover this decrease, the stronger J' is needed 
for the AF LRO to appear, yielding an initial increase of 
J' c {x). If, however, our scenario is applied to the system 
just at J' = J' c , where the correlation lengths ^ v are 
infinite, it is expected that the AF LRO is induced by 
an infinitesimal concentration of dilution, since {J mn } 
is quasi-long-ranged even for a vanishingly small x and 
the prefactor in eq. (13), jTo, depends more weakly on J' 
than J a f (oc J' 2 ) does. This means that J' c {x) has a neg- 
ative slope at least at x — + . Actually, from eqs. (12) 
and (13) combined with the critical behavior of £^' y at 
J' — Jc pi the decrease of J' c (x) at x ~ is deduced. 
However, quantitative details of its behavior in the weak 
dilution limit arc beyond a scope of the present work. 

5. Conclusion and Discussion 

In the present paper, we have investigated the bond- 
dilution effects on the non-magnetic ground state of the 
2D quantum AF Heisenberg model consisting of bond- 
alternating chains by means of the QMC simulations, and 
have proposed a scenario that the quantum phase tran- 
sition induced by bond dilution is originated from the 
competition of the two effective AF interactions {J mn } 
and J a f. In particular, the proposed simple mapping of 
a bond-diluted system to a non-diluted system catches 
up the essential mechanism of the bond-dilution-induccd 
AF LRO in the system so long as the effective interac- 
tion J(x), which represents an averaged value of { J mn }, 
is less than or comparable to J a f, or in other words, J' 
is not too small. 

In our arguments so far described, we have kept away 
from an interesting problem whether the disordered state 
in the present system may involve a novel aspect of quan- 
tum random systems such as the quantum Griffiths (QG) 
phase or not. By our scenario with the averaged inter- 
action J(x) we completely neglect this possibility. For 
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J 1 = 0, i.e., the system is exactly one-dimensional, the 
system with < x < 1 is known to be in the QG phase 
characterized by finite correlation lengths but with a con- 
tinuous distribution of spin gaps up to zero. 31-34 If this 
QG state is unstable against the AF LRO state due to 
the introduction of J' even of an infinitesimal magnitude, 
the latter state invades to x — > 0, i.e., a reentrant tran- 
sition from the disordered state to the AF LRO state as 
decreasing J' is expected for x smaller than 0.11. Also for 
a moderate value of J' for which we have observed the 
quantum phase transition from the disordered state to 
the AF LRO state, there remains the problem of a pos- 
sible existence of the QG state. Actually in our prelim- 
inary analysis we observed in the J' = 0.3 system that 
A^ 2 ) of cq. (2) decreases much faster than the inverse 
of the spatial correlation length as x increases. 17 This 
strongly suggests an existence of the QG state. However, 
in order to confirm a continuous distribution of spin gaps 
up to zero, we have to carry out QMC simulations at 
correspondingly low temperatures. This problem as well 
as that of the possible reentrant transition mentioned 
above remain as challenging future works in the field of 
the computational physics. 

In order to distinguish experimentally the charac- 
teristic difference between the site- and bond-dilution- 
induced AF LRO states pointed out in the present work, 
one needs a suitable material; a quasi-one-dimensional 
Heisenberg antiferromagnet with small enough J' and 
rigid bond alternation, and without frustration. Fur- 
thermore these properties are desired not to change 
significantly by dilution processes. In the spin-Peierls 
compound CuGei-zSizOa, which is effectively a bond- 
disorder Heisenberg antiferromagnet, 2 it was experimen- 
tally observed that the effective spins are induced not 
near Si impurity sites but rather at sites between them 
in the AF LRO phase. 9 The occurrence of the AF LRO 
is then attributed to the nonlocal reconstruction of bond 
alternation by Si substitution. To understand this type 
of phase transitions in the spin-Peierls state, we have to 
analyze models with proper spin-lattice coupling, which 
will be another challenging numerical work in near fu- 
ture. 
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